submitted to J. Stat. Phys. 

Escape-rate formalism, decay to steady states, 
and divergences in the entropy-production rate 

Jiirgen Vollmer, 1 '!^] Laszlo Matyas, 2 and Tamas Tel 2 '[] 

1 Max- Planck- Institute for Polymer Research, 
Ackermannweg 10, 55128 Mainz, Germany 
2 Institute for Theoretical Physics, Eotvos University, 
P. O. Box 32, H-1518 Budapest, Hungary 
(Dated: December 14, 2001) 



In summer 1997 we were sitting with Bob Dorfman and a few other 
friends interested in chaotic systems and transport theory on a terrace 
close to Oktogon in Budapest. While taking our (decaf) coffee after a 
very nice Italian meal, we discussed about logarithmic divergences in 
the entropy production of systems with absorbing boundary conditions 
and their consequences for the escape-rate formalism. It was guessed 
at that time that the problem could be resolved by a careful discussion 
of the physical content of the absorbing boundary conditions. To our 
knowledge a thorough analysis of this long-standing question is still 
missing. We dedicate it hereby to Bob on occasion of his 65th birthday. 
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I. INTRODUCTION 

The escape-rate formalism [Q, |^, [| aims at identifying transport co- 
efficients based on the asymptotic decay rate of an initial non-stationary 
density profile towards an empty steady state selected by absorbing bound- 
ary conditions where all particles disappear from the interior of the system. 
In this formulation the relaxation problem has widely been studied in the 
context of Markov chains (c/. the sections on survival probabilities in 
and on absorbing states in j6|), as well as for deterministic chaotic systems 
(7[ ||, [|. However, the choice of an empty asymptotic state places severe 
constrains on the formalism. In particular, the thermodynamic entropy 
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production picks up logarithmically diverging contributions at the bound- 
aries. After all, a vanishing density is physically unrealistic. In particular, 
it leads to the breakdown of the concept of an entropy-production density 
even in the framework of classical irreversible thermodynamics. Close to 
a boundary at x = the absorbing boundary condition requires a density 
profile of the form p(x) — ax. For purely diffusive particle transport, a 
neighborhood of size A gives then rise to the entropy production 1 



(A) = /° +A dx D ( ^B) p = Da km In $. 
Jo V P J s^o 5 



(1) 



In the present paper this logarithmic divergence will be discussed from 
the point of view of spatially extended chaotic systems whose transport 
properties fully agree with the predictions of irreversible thermodynamics. 
To keep the calculations as transparent as possible, the discussion is given 
for isothermal multibaker maps. However, the described picture should 
apply in general. 

In Sect. H Eq. (|l|) is contrasted with the prediction of the irreversible 
entropy changes in an entropy based on the conditional density character- 
izing the chaotic saddle forming the backbone of transport in the system. 
This prediction always yields finite values. The isothermal multibaker map 
is introduced in Sect. Ill, where also its entropy balance is worked out. 
Sect. [V deals with the normal modes of the coarse-grained time evolution. 
This allows us to address in Sect, [v] the origin of the logarithmic divergences 
in the entropy production from the point of view of a microscopic reversible 
dynamics. The analysis makes use of the eigenvalues of the time-evolution 
operator ^, [h]], which does not depend on the nature of the asymptotic 
state. To underline this observation, we also discuss asymptotic states with 
uniform nonzero densities, and point out how the divergences are lifted by 
an arbitrarily small background density. The discussion also shows why the 
changes of the irreversible entropy based on the conditional density differ 
from the thermodynamic expectation. The presence of an arbitrarily small 
background density in the asymptotic case turns the entropy production to 
be finite, but different from the prediction based on the conditional density. 



In the concluding Sect. VI these findings are complemented by a discussion 
of the behavior of the entropy production in systems relaxing towards a 
typical nonempty steady state of finite density as compared to the case of 
small (or even vanishing) background densities addressed in the main part. 



1 The Boltzmann constant is taken to be unity throughout this paper. 
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II. ENTROPY PRODUCTION BASED ON CONDITIONAL 
INVARIANT MEASURES REVISITED 



One of the early studies of thermodynamic entropy production in deter- 
ministic dynamical systems was based on the escape-rate formalism intro- 
duced by Gaspard and Nicolis As a generalization of it, Breymann, Tel 
and Vollmer Q considered open dissipative dynamical systems in continu- 
ous time. To characterize their irreversible features, they suggested to use 
the entropy 

a(i) = - [dx V (t) 0) In ifj {t) (x) (2) 



based on the normalized conditional phase-space density tp ^'(x), describ- 
ing the probability to find a point which has not yet escaped the system 
by time t at phase-space coordinate x. Because tp^ (x) is a single-particle 
property, s(t) can be considered as a specific entropy (total entropy per 
number of particles). The initial condition is selected from an arbitrary 
smooth distribution ip^(x). As time goes on, the phase-space volume 
of ipW is exponentially shrinking as exp [— q,(x) t], where the local phase- 
space contraction rate <;(x) = \i(x) — n is a smooth function of the 
coordinates 0. Here, k is the escape rate from the system, and Xi(x) 
are the local Lyapunov exponents that describe the phase-space contrac- 
tion in the independent directions i in phase space (c/. |?], || for detailed 
discussions of conditional invariant densities and Lyapunov exponents). 
On its support the value of the conditional density ip^(x) is exponen- 
tially increasing due to its normalization. More precisely, it increases like 
npit+dt) (x) — exp [s(x) dt] ip^(x) x^^Hx), where x ( - t+d *^ is the character- 
istic function of the support at time t + dt. The entropy s(t + dt) at time 
t + dt can be determined by inserting this relation into Eq. (Q) : 

s(t + dt) = -dt [ dxs(x)il>®(x) e< x ^ dt x it+dt) (x) 



dx iP (t) (x) In U^(x)x it+dt) (x)\ e^ x)dt x {t+dt) {x). (3) 

In both integrals the decrease of the support of tp is counterbalanced by 
the factor exp[q(x)di\. The first integral is the phase-space average <f of 
s(x). For dt — > the second one tends to the specific entropy s(t) at time 
t. Hence, in the long-time limit the time derivative of the entropy 

S=-^5>- k ' (4) 

i 

is the difference of the sum of the average Lyapunov exponents on the 
saddle, and the escape rate k from the saddle. The average is taken with 



4 



J. Vollmer, L. Matyas, and T. Tel 



respect to the density ip(x) of the conditionally-invariant measure. This 
measure is time independent. Its support is the unstable manifold of the 
chaotic saddle, i.e., the union of the never escaping orbits in the system. 
The fact that the time derivative approaches a constant reflects the ever 
refining fractal structures in the density due to the chaoticity of the dy- 
namics. 

We now compare s(t) with a coarse-grained entropy s^° s \t) computed in 
an analogous way from a coarse-grained conditional density ^^'(t), which 
— in contrast to ip — does converge towards a stationary distribution. The 
irreversible entropy production is then obtained as (cf. 

p(irr) = d j s (cg)^ _ s ( t )j l-£i5- K-^Xi. (5) 

i 

It measures the lack of information on the microscopic state due to the 
finite resolution of the coarse-grained description. Similarly to s, p( ur ) is 
a specific quantity. 

In systems with a reversible dynamics the phase-space contraction is 
proportional to the displacement in the direction of an applied field |n], [l2j 
[13[ [14| |l5| , |l6f . Consequently, in an open system with reversible dynamics 
the sum of the average Lyapunov exponents on the chaotic saddle is zero 
since the average number of steps in the direction of positive phase-space 
contraction is the same as in the opposite. Therefore, its specific irreversible 
entropy production amounts to the escape rate, 

= k. (6) 

In the following we revisit this argument in the light of recent devel- 
opments |l3| , |l7| dealing with steady states instead of empty asymptotic 
states. We work out the irreversible entropy production for an isothermal 
multibaker map with reversible microscopic dynamics subjected to absorb- 
ing boundary conditions. 



III. THE ISOTHERMAL MULTIBAKER MAP 



Multibaker maps model particle transport in spatially extended systems 
by a chain of mutually interrelated baker maps (l^, [l4|, [l7| [l^, 

ST 



24|. They consist of N identical cells of width a and height 1 



(the phase-space) in the (x,p) plane. The cells are labeled by the index m 
(Fig. la). After each time unit r, every cell is divided into three columns 
(Fig. lb). Here we consider the case when the right (left) column of width 
ar (al) is mapped onto a strip of width a and of height I (r) in the right 
(left) neighboring cell. The middle one, which is of width as, preserves its 
area, such that its image attains a height s, and r + I + s = 1. There are 



Escape-rate formalism, steady states, and entropy production 



5 



(a) 



CO 
TD 
C 

o 
-Q 



I , im i i i rren 



m = 



a 
1 



a 
2 



a 
N 



- CD 2 

c -a 

_ , =3 C 

o o 

La u 

N+1 



(b) 



(Pm-l) 



(Pm+1> 



1 ar 1 1 as 1 'al' 
m-1 I m I m+1 



a 
m 



!JPm,r) 
(Pm'.s ) 

(Pm,l ) 



FIG. 1: Graphical illustration of the action of the multibaker map of length 
L = aN on the phase space (x,p) over a time unit r. (a) The mapping is defined 
on a domain of -/V identical rectangular cells of size a x 1, with boundary condition 
imposed in two additional cells and N + 1. (b) The action of the map in any 
of the cells over time unit r is illustrated by the deformation of the labels R, S 
and L in the three branches of the map ending up in cell m. The average value 
of the density on the cells (strips) [cf. Eq. ([?])] is given on the margins. 



more general parameter settings conceivable, but earlier work ||l3| , |l6| , |l7| 
showed that the associated macroscopic behavior is then not compatible 
with irreversible thermodynamics. 

The dynamics of the multibaker map models a microscopic dynamics 
described in the single-particle phase space. It is deterministic, invertiblc, 
chaotic, and mixing [ p~5 , p5[ . To describe irreversible processes one follows 
the coarse-grained densities p m obtained by averaging over the cells Jl3| , 
~ 3, 17 . To emphasize the particular choice of coarse-graining over the 
Is, the coarse-grained densities are also called the cell densities. The 
dynamics of the multibaker map is the same for all cells. There might be 
inhomogeneities in the densities, but the evolution equations are translation 
invariant. 



E 

ce 



A. Evolution of the cell density 

In order to find results consistent with non-equilibrium thermodynamics 
we always consider initial conditions with a uniform density in every cell 
to. This is convenient from a technical point of view, and does not lead to 
a principal restriction of the domain of validity of the model as discussed 
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in |b|, [16| ^6|. Under such conditions the parameters r and I can be con- 
sidered as transition probabilities from a cell to its right and left neighbor, 
respectively. After one step of iteration the densities p' m i on the strips 
i = B,, S, L of cell m are (c/. Fig. |l|b) 



r „ J J 1 



Pm,r ; Pm—li Pm.s — Pmt Prn,l ~ Pm+l- (7) 

The factors r/l and l/r give rise to local contraction or expansion of the 
phase-space volume. One of the factors is larger than unity, and character- 
izes a local contraction, while the other gives rise to an expansion. 

Since in every cell the density remains uniform in the horizontal direc- 
tion, this update holds at all times such that the coarse-grained density 
distributions p m and p' m = rp m -\ + sp m + lp m +i at the respective times n 
and n + 1 are related by the master equation 

Pm = (1 - r - l)p m + rp m _i + lp m +l. (8) 

Multiplying the equation by r _1 and introducing the current 

jm = -(rp m ~ I Pm+l) (9) 
T 

through the right boundary of cell to, Eq. (g) appears in the form of the 
continuity equation 

Pm ~ Prn jm ~ Jm — 1 (10) 

t a 

The current through the left boundary of cell to is the same as the current 
flowing through the right boundary of cell m — 1. 



B. Diffusion, drift, and the macroscopic limit 

The transition probabilities r and I govern the evolution of the coarse- 
grained density p m . In view of the master equation (j^), the cell-to-cell 
dynamics of the model is equivalent to the dynamics of an ensemble of 
random walkers with fixed step length a and local transition probabilities r 
and I over time unit r. In terms of the local drift v and diffusion coefficient 
D p^] the transition probabilities r and I can be expressed as 




such that the current appears in a form very close to its thermodynamic 
counterpart, viz. 



(12) 
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The macroscopic limit expresses a separation of scales where density 
gradients inside cells may be neglected, while density differences between 
the cells and the temporal evolution of the cell densities are only taken 
into account in leading order. In other words, in this limit a <C L = aN , 
t <C L 2 /D, L/v, and upon introducing the quasi-continuous spatial and 
temporal variables x = am and t = rn, the current j m [Eq. ( |12|) takes the 
macroscopic limit j(x;t) = v p(x;t) — Dd x p(x,t) , while Eq. (|1C) becomes 
dtp(x; t) = —d x j(x; t) = —v d x p(x; t) + D d%p(x; t) . The macroscopic den- 
sity evolves according to the advection-diffusion equation. 



C. Local entropy balance 



The coarse-grained entropy of cell m is defined as 

S m = -a/9 TO ln-^-. (13) 
P 

It fulfills a local entropy balance in direct analogy to the one of irreversible 
thermodynamics fl3, 14 1. In this equation p* is a constant reference density 



that is introduced for dimensional reasons. In classical physics it expresses 
the free choice of the origin of the entropy scale. 

To derive the balance equation for (Il3|) one identifies at any given time 



(G) — (G) 

the difference S m — Sin of the coarse-grained and the Gibbs entropy Sm 
as the information on the microscopic state of the system which cannot 
be resolved in the coarse-grained description. The Gibbs entropy is an 
analogous expression to ([l3]) given in terms of the non-coarse-grained phase- 
space density p, as —phip/p* integrated over the cell. For a coarse-grained 

(G) 

initial distributions the entropies coincide initially (i.e., S m — Sm )■ After 
one time step the entropies become (c/. Fig. |l]b) 



c(G)' 



sp m In — + rpm-i m — — — + lp m +i hi -- 



p* \l p* J \r p 



(14) 



and 



^ = -ap' m ln^. (15) 



The temporal change of the lack of information is identified with the 

irreversible entropy production AjiS m , and the change (Sm — Sin ) of the 
Gibbs entropy with the entropy flux A e S m - Thus, one obtains the discrete 
entropy balance in any cell 

Sm Sm ^eS m , ^iS m r-\n\ 

= + . (16) 
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The form of the entropy production is [cf. ( |l4| ) and (15)] 



T 

a 

T 



.-(G) 1 



-p' m In — + rp m -x In ( '- 



Pn 



r Pm-l 



I p n 



, I / Pm+1 

lpm+i In 

Prn 



,(17) 



and the entropy flux becomes 
A c c(G)' _ q(G) 

"-^e^m ^>m ^>m 



{p'm ~ Pm) m ~T + r /°m-i ln 



r p m -i 
I Pm 



I Pm+1 In 



r Pm /. 



J) 



In the macroscopic limit all expressions reduce to the respective predictions 
of non-equilibrium thermodynamics jl3|, Note that the entropy 

production does not depend on the choice of the reference density p* . The 
density of irreversible entropy production is then AiS m /(aT). 



D. Global entropy production 

The escape-rate formalism addresses the balance of the global entropy 
of the chain. The global coarse-grained entropy is 

N N 

Stot = ^ S m = p m \n-^. (19) 

m— 1 m—1 

The associated global entropy production rate along the chain is £( 1IT ) = 
Sm=i AiS'm/r, and the total specific irreversible entropy production is 
obtained as 

y(irr) AC 
p(irr) _ _ l^m=l , on \ 

N ~ tN [ } 

where TV = Y] —1 ap m is the number of particles in the chain at time nr. 
The total entropy production can be rearranged to take the form 

V(irr) a ' , Pm , a ( , r P0 , , ^PjV \ 

S v ; = > p m ln — + - [rpo In IpN+i In -. 

T m~^l P m T V 'Pi 'PW+l / 

N-l 

H — /j(TPro - lp m +i)ln- — — , (21) 



m—1 
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where p an d Pn+i are the densities in the boundary cells. It can conve- 
niently be split into four terms £( irr ) = sf rr) + s[ irr) + S^ r) + E^. The 
first one 

T n Prn 

is the contribution from the temporal evolution of the density. The con- 
tribution proportional to ln(r/Z) of the last term in ( ^T| ) contains the irre- 
versible entropy production 

= M- In - - M— (23) 

r / L> 

due to the presence of the drift v. It does not depend on the particular 
density distribution so that we could immediately specify its macroscopic 
limit (indicated by — >). By means of (§) and (|g) the remaining part can 
be written as a sum of two terms. One of them, 

N 
m— 1 

characterizes the contribution of mixing of the neighboring densities. In 
order to arrive at this form the ratio p m /pm+i of the densities appearing 
at the right hand side of ( |2l| ) was written as [(p m / P*)/ {Pm+i/ P*)]- The 
rest 



v (irr) _ a 



, r Po . , , IpN+i , , Ipi , rp N 

rp Q In — + /pAr+i In Ipi In — - - rp N In — - 

tp* rp* rp* (p* 



(25) 



yields the boundary contribution. We shall be interested in the difference 
between £( irr ) and S^ lrr \ called the irreversible entropy production S^ciaL 
connected to the relaxation process, 

E (xrr) ^ E (irr) ^(irr) ^(irr) ^ 
relax t 1 "mix 1 b \^"J 



IV. NORMAL MODES OF THE COARSE-GRAINED TIME 

EVOLUTION 

A. Decaying modes and the steady state 

(n) 

We arc interested in the evolution of the density distributions p m sub- 
jected to a fixed constant boundary condition 

Pq - Pn+i - Pb (27) 
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at any time step n. Asymptotically, p^m always approaches the uniform 
density pb- The time evolution of the density can be explored by expanding 
the deviation — pg of the initial distribution from the asymptotic state 



in terms of normal modes 



>](r 



They vanish at both boundaries, and 



only change in amplitude but not in their shape, 



rM(n+l) 
> m 



exp(-7 I/ r). 



(28) 



The integer v labels different modes. There are as many independent modes 
as the number N of the cells v = 1, . . . , JV. The normal modes take the 
respective forms 



exp 



(-7„nr) (j 



%/2 



N + 1 



(29) 



Substituting the ansatz into Eq. (g) and rearranging the trigonometric 
terms, one finds the decay rates 



7k = --ha 

T 



1 - (r + I) + 2Vrl cos 



N + l 



7T 2 D 



m < 30 > 



For a general initial condition the asymptotic decay is governed by the 
slowest non- vanishing decay rate, 71. It coincides with the escape rate k of 
the transiently chaotic motion || inside the chain (i.e., k = ji) |l|, ^8| . 

The macroscopic limit of the decay rates has a clear physical content. For 
vanishing v it states that relaxation is related to the typical diffusive decay 
rate D/L 2 of structures of size L. The factor ir 2 characterizes the geometry 
of the considered region (a band of width L with straight, parallel walls in 
our example). More complicated geometries have been studied recently by 
Gaspard |18 , and Kaufmann and collaborators 123, BQ]. 



For a biased motion v ^ the drift singles out one side of the system 
and sweeps out the particles in that direction. This mechanism dominates 
when the time L/v to cross the system by the biased motion becomes 
shorter than the typical time scale L 2 /D of diffusion, i.e., for 



Pe: 



vL 
~D 



(31) 



much larger than unity. In the context of hydrodynamics, Pe is called 
Peclet number It measures the importance of diffusion relative to 

advection. Strong diffusive effects are indicated by small Peclet numbers. 
For fixed finite v and D, the Peclet number is always large for a sufficiently 
large system size L. 
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B. Long-time relaxation and the slowest mode 

For sufficiently long times n ^> 1, the coarse-grained density closely 
approaches the first normal mode. Therefore, the density can be expressed 
as 

^ PB + (N^ - N^) ^n, (32) 

where A/"*- 00 - 1 = pbL is the particle number in the background which is also 
the asymptotic particle number in the system, and 

_4 ,nm/2-(JV+l)/4 m7T 

^ = lllJ sin (33) 

is the coarse-grained conditionally-invariant density. It is normalized to 
unity (1 = a X!m=i Vw)> by virtue of the normalization constant A, which 
is invariant under the exchange of r and I. Carrying out the summation of 
the complex geometric series defined by (|33| ) one finds 

1 — expf— kt) 1 1 

A = L ~ 



(f) (Ar +1 )/4 +(f) -(Ar +1) /4 sin _^_ 

2 cosh^ 27rDcosh^' V ; 

Here, the relation (r/l) N ^ 4 = (l + av/D) N ^ 4 — > exp(Pe/4) has been used to 
evaluate the macroscopic limit. The asymptotically decaying density takes 
then the form 

p(x,t) =PB + j (A^> - exp (Pe^^) sin^. (35) 



V. BOUNDARY CONTRIBUTIONS TO THE IRREVERSIBLE 
ENTROPY PRODUCTION 

A. Absorbing boundaries 

In the case of a long-term relaxation towards an empty state (po = 
Pn+i = Pb = 0), Eq. ( pq ) holds for the full density p m , and one can write 
[see Eq. @] 

S t (irr) ee nj\f exp(-KT) = kAT' -> kJV, (36) 
where N 1 is the number of particles at time (n + 1)t. 
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The mixing term (E4T) can be expressed by means of the total entropy 



(19) to obtain 



^—(e— -1)5* 



— « St. 



(37) 



and in view of po = Pn+i = Pb = the boundary contribution (p5| 
becomes 



,(irr) 



a ^4 „• 7r 



,(iV-l)/4 



hi 



AN /r\(^+3)/4 . w 



+'(?) 



-(JV-l)/4 



hi 



£(*)' 



JV+1 



(JV-l)/4 



+ *(?) 



-(JV-l)/4 



hi 



44 sin ^ 



r ( r ) t^" 1 )/ 4 - I (r ) "^"^l iY±3 lQg r . (38) 



Observing that in the macroscopic limit both r and I are in leading order 
equal to rD/a 2 , and that log(r//) — > av/D, one obtains for the relaxation 
contribution to the total irreversible entropy production 



,(irr) 
J relax 



Pe , Pe , 

1 tanh In 

4 4 



/ AT n 2 1 + (Pe/27r) 2 a 
\Tp*Y cosh(Pe/4) L 



reSt, 



(39) 

This expression shows the expected logarithmic divergence since (a/ L) — > 
in the macroscopic limit. On the other hand, the result can properly be 
interpreted only after evaluating Stot- In particular, the reference density 
p* has to drop out again in the final result. 

In the macroscopic limit the sum over m in the definition of Stot becomes 
an integral. By using Eqs. (|32]) one finds 



-AN j dx exp (-Pe 
Jo 



2x-l 



In 



AN 
Lp* 



exp 



i j sin(7r:r) 

2x-l 
-Pe 



sin(7rx) 



(40) 



Applying the relation (|34), we see that the specific entropy 

S to t=Nf(pe,^) (41) 

is a function of the Peclet number and of the ratio of the average density 
in the system N/L and the reference density p*, i.e., of another dimen- 
sionless constant that involves the parameter p* selecting the origin of the 
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entropy scale (which, as mentioned earlier, is an arbitrary number in clas- 
sical physics). The total entropy cannot be evaluated exactly. However, to 
obtain its behavior in leading order for very large and small Peclet num- 
bers it is sufficient to approximate the expression under the logarithm by 
its maximum value. In the two limiting cases one thus finds in leading 
order in M 

{ -Nm ( Pe^-) for Pe > 1, 
Stot = { r ) (42) 
1 -AHnf^y for Pc< 1. 

This implies for the specific entropy production 

r > E (i 7 ) fKlnfPe" 1 ^) for Pe > 1, , 

p(iTi) _ relax _ J V a I ^ ' (43) 

TV \ Khi | for Pe< 1. 

The result shows the expected logarithmic divergence due to the boundary 
terms. It should be considered as an example clearly showing the inade- 
quateness of global quantities for characterizing thermodynamic properties 
[cf. Eq. (Q)]. The reason for the breakdown of the prediction (||) lies in 
the fact that the argument leading to this result focuses on the shrinking 
of the support of the measure by assuming the smoothness of the distri- 
bution along the unstable manifold. It thus entirely disregards that the 
density is very inhomogeneously distributed as a consequence of the ab- 
sorbing boundaries. It should also be noted that the result obtained for 
the diffusive case Pe <C 1 is the analog of the thermodynamic expression 
since in this case p = Nir/(2L) sm(xir/L), such that the parameter a 
of Eq. (0) takes the value a = Afn/(2L 2 ) and Da = Nn/2. The presence 
of factor 1/2 is due to the fact that Eq. gives the contribution of one 
end only. 



B. Influence of a small background density 

We now assume that ps is nonzero but much smaller than p m except for 
a narrow boundary layer where the sine of ([33|) approaches zero. In that 
case Eq. (|8|) gives an upper bound to the ratio of densities at successive 
times, that is very accurate in the interior of the system. Consequently, the 



evaluation of (26) carries over except that the boundary term s[ lrlJ picks 



up contributions due to the finite density ps in the cells and N + 1. In 
the macroscopic limit this term becomes 

S[ irr) -f -k{M - *M ) ( 1 + In BE) « -kN ( 1 + In BE. ). (44) 
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Since the total entropy Stot & n d the contribution sj 1 "' 1 to the entropy 
production do not significantly change for a sufficiently small background 
density, the full entropy production becomes in the macroscopic limit 

^l = -^ln^-^ tot . (45) 



After substituting 5 tot from (42), we obtain the specific irreversible entropy 
production for Af > Af^ , 

, v f KlnfPe^-H^-r) for Pe > 1, 

p(irr) = I V JV(°°> ) , 4 gx 

\ Kln(^y) for Pe< 1. 

The result clearly shows that the logarithmic divergences in the entropy 
production of the previous case are due to the vanishing of a physically 
indispensable background density ps- 



VI. DISCUSSION 

The result (^6|) involves only well-behaved macroscopic quantities, and 
the logarithm of the ratio of the number Af of particles in the system 
over the number TV*- 00 -* = psL of particles approached in the steady state. 
At intermediate times, where Af is still much larger than Af^°°^ , the ratio 
Af /Af(°°) decreases to a good approximation exponentially like exp(— nt) 
such that the rate of irreversible entropy production starts to decrease lin- 
early like —K 2 t. During this time regime the boundary contribution is by 
a factor of log(Af /Af^) larger than the bulk contributions accounted for 
by Eq. (|J). Hence, even in the more realistic setting accounting for a finite 
background density, Eq. (|J) only describes a sub-dominant contribution to 
the entropy-production rate. The reason for its failure is that the contribu- 
tions arising from the spatial distribution of the particles and the induced 
inhomogeneous particle currents are not adequately taken into account by 
an entropy based solely on the conditionally-invariant measure. 

In spite of the strong contributions due to the boundary terms, how- 
ever, the entropy production remains proportional to the escape rate k 
even in more general situations. Besides for the short times, where this 
condition follows from Eq. (f46|), this can be easily illustrated in the long- 
time regime, for which (Af — Af^^ijj = AAftp <C pb- Consequently, the 
entropy-production rate is obtained as 



D V P(a) 
— f dx p(x) — — D [ dx ip(x) d^.ijj(x) 

D J PB J 



Escape-rate formalism, steady states, and entropy production 
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*X* +K Wr g <Pe) (47) 

where <?(Pe) — J dx ip 2 (%) is a function of Pe only. The corresponding irre- 
versible entropy production due to relaxation is £|."iax = £<( lrr ) — N ' v 2 / D. 
In the second line in Eq. ( |47| ) the term proportional to d x tp does not appear 
since its integral vanishes. Moreover, an integration by parts was used to 
obtain a second spatial derivative of the density which, according to the 
advection-diffusion equation, is proportional to its time-derivative, i.e., it 
amounts to — nil) for the slowest decaying mode ip (again it is used here 
that terms proportional to d x ^p and ijjd x ^j vanish under the integral). The 
function g(Pe) can easily be evaluated for the multibaker map, but in gen- 
eral it depends on the shape of the system. Thus, the specific irreversible 
entropy production taken with respect to AJ\f is 

P^ = ^=«j^9(Fe). (48) 



Equation (48) implies that even in a general thermodynamic setting the 
relaxational entropy-production is proportional to the escape rate k, which 
characterizes the approach towards the stationary state. In contrast to the 
dynamical-system arguments based on the escape-rate formalism, the 
term involves in general a non-trivial function g(Pe) of the Peclet number, 
and it has an amplitude AjV /A/"*- 00 ) that is exponentially decaying like 
exp(— nt). 

It will certainly be interesting to investigate more closely the connection 
between the escape-rate formalism and the decay to systems supporting 
non-trivial stationary states. Another first step in this direction, which 
complements the present approach was suggested by Gilbert et al |32| ], 
who recently discussed the approach towards equilibrium in a system with 
periodic boundary conditions. 
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